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ABSTRACT 


Recent events in the United States of America and Pakistan have exposed the 
shortcomings of existing planning in relief and humanitarian assistance in the face of 
large-scale natural disasters. This thesis develops a two-stage stochastic optimization 
model to provide guidance in the pre-positioning of relief units and assets, where budget, 
physical limitations and logistics are taken into account. Stochastic data include the 
number of survivors in each potential affected area (AA), the amount of commodities that 
needs to be delivered to each AA and the transportation time from each relief location 
(which reflects scenarios where, for example, roads are blocked). As first-stage decisions, 
we consider the expansion of warehouses, medical facilities and their health care 
personnel, as well as ramp space to facilitate aircraft supply of commodities to the AAs. 
The second-stage is a logistic problem represented as a network, where maximizing 
expected rescued survivors and delivery of required commodities are the driving goals. 
This is accomplished through land, air and sea transportation means (e.g., CH-53 
helicopters configured for rescue missions), as well as relief workers. The model has 
been successfully assessed on notional scenarios and is expected to be tested on realistic 


cases by personnel who are involved in relief planning. 
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EXECUTIVE SUMMARY 


Recent events in the United States of America and Pakistan have exposed the 
shortcomings of existing planning in relief and humanitarian assistance in the face of 
large-scale natural disasters. Shortcomings include factors like poorly pre-positioned 
relief units and assets, insufficient budget allocated, and poor post-disaster relief 
planning. These contributed to long response time, which in turn caused the delay of 
medical evacuation and supply of commodities, and perhaps the death of hundreds of 


potential survivors. 


We develop a two-stage, linear, mixed-integer program called Pre-positioning 
Optimization Model (POM) to provide guidance in the pre-positioning of relief units and 
assets, where budget, physical limitations and logistics are taken into account. Stochastic 
data include the number of survivors in each potential affected area (AA), the amount of 
commodities that needs to be delivered to each AA, and the transportation time from each 


relief location (which reflects scenarios where, for example, roads are blocked). 


As first-stage decisions, we consider the expansion of warehouses, medical 
facilities and their health care personnel, as well as ramp space to facilitate the supply of 
commodities by aircraft to the AAs. These decisions are termed first-stage as their 
implementation must be carried out well in advance, before a disaster strikes. The 
second-stage is a logistic problem represented as a network, where maximizing expected 
rescued survivors and delivery of required commodities are the driving goals. This is 
accomplished through land, air and sea transportation means (e.g., CH-53 helicopters 
configured for rescue missions), as well as relief workers. An unmet commodity penalty 


is built in the objective function to penalize the total number of rescued survivors. 


We implement POM in the General Algebraic Modeling System and use CPLEX 
as the solving engine. Two hypothetical test cases (Earthquake and Hurricane) are used 
for testing. To simulate the different levels of disaster severity, each test case has its 
distinct set of scenario data (e.g., number of potential survivors, demand of commodities 


and impassable roads) and a hypothetical probability. A bound on the optimal stochastic 


XVil 


solution is obtained by solving POM with perfect information (i.e., as a weighted average 
of the deterministic solutions for each scenario). The optimal stochastic solution is also 


compared with a heuristic solution based on the optimal plan for an “‘average scenario.” 


The benefits of solving POM using stochastic optimization are apparent: total 
number of rescued survivors (after penalties for unmet commodities are accounted for) is 
improved by 11% and 37% in the earthquake and hurricane test cases, respectively, with 
respect to the heuristic solution. But even if a deterministic scenario were considered, 


finding an optimal solution without POM would be nearly impossible. 


POM can be used as a tool to aid in planning budget allocation, sitting of relief 
forces, and pre-positioning of warehouses and other assets. We recommend further study 
using validated data. Our tool is expected to be used by personnel who are involved in 


relief planning. 
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I. INTRODUCTION 


The term disaster usually evokes images of massive material damage and great 
human distress caused by some swift catastrophe. A disaster can be defined as a sudden 
event that disrupts the social structure, and prevents execution of some or all of the social 
structure’s essential functions [Foster, 1983]. Disaster may surprise rich and poor alike, 
and show how even highly capable governments are unprepared to deal with the massive 
exigencies of emergency relief. Recent natural disasters in the Pakistan [Wikipedia, 
2006] and the United States of America [Wikipedia, 2005] exposed the shortcomings of 
existing planning at both domestic and international cooperation levels. In the above 
examples, the affected area’s authorities lacked the manpower and logistic preparation to 


assemble and coordinate sufficient teams to aid in the humanitarian rescue. 


Humanitarian planning and assistance is a large-scale process that requires the 
collaboration and coordination of government agencies and other organizations. For 
example, in the U.S., participants include the Federal Emergency Management Agency, 
the National Guard and other military forces [U.S. Army, 2005], as well as the American 
Red Cross and Food for the Hungry [OPM, 2006]. These teams often operate in a 
complicated environment that includes temporally and spatially varying demand patterns, 
administrative, legal, and political constraints; and ill-defined objectives [Chaiken and 
Larson, 1972]. Unpreparedness results in the inability to handle “chaos,” which in turn 


increases the loss of lives. 


In an emergency response situation, it is critical to minimize the response time, 
because each minute wasted diminishes the chances of rescuing potential survivors. 
Therefore, strategic positioning of relief units (e.g., medical facilities and health care 
personnel) and assets (e.g., CH-53 helicopters configured for rescue missions) becomes 
crucial. Unfortunately, it is impossible to place relief units and assets (RUAs) 
everywhere. A new 464-bed hospital costs around $550 million (Department of 
Medicine, John Stroger Hospital, 2004), a general doctor costs around $420 per day 
[Yam, 2004], a new CH-53 cost $25 million in 1993 (FAS 1999), and a new 30,000- 


square-foot (approximately half the size of an American Football pitch), one-story 


warehouse costs around $2.3 million [Department of Commerce and Economic 
Opportunity, 2004]. It is clear that we can afford only a limited number of each of these 
assets. Common sense dictates to position RUAs “close to the probable affected areas” 
where a natural disaster may strike. On the other hand, it also appears reasonable to 
scatter the RUA locations to enable partial coverage of multiple affected areas by a RUA 
location, and simultaneous coverage of each possible affected area by multiple RUAs 


(e.g., to avoid dependency on a single RUA in case it is not accessible after the disaster). 


During the first few days after a disaster, influx of logistics also poses a challenge. 
Logistic infrastructure requires pre-positioning of warehouses and ramp space for aircraft 
to deliver massive amounts of commodities to the affected areas. “Emergency/Relief 
Logistics is the basic task of a logistics system: to deliver the appropriate supplies, in 
good condition, in the quantities required, and at the places and the time they are needed. 
Relief Logistics encompass the relocation of disaster affected people, transfer of 
casualties, and the movement of relief workers” [Hanaoka and Qadir, 2005]. These 
logistical implications have been studied by many researchers, such as Darcy [2005], who 
discusses the problems in accessing the affected areas in the Indian Ocean tsunami crisis 


in 2004. 


Modeling the needs before and after a disaster is a complicated process and the 
complication becomes increasingly intricate by the large amount of RUAs involved. For 
example, the transportation means can be land, air and/or sea. The number of 
transportation means needed will hinge on the budget. Expansion of warehouses, medical 
facilities, ramp space for aircraft and many more considerations are factors that will 
affect the budget allocation directly. This thesis develops mathematical models to help 


strike a balance among these considerations. 


The analysis focuses on the strategic level planning of humanitarian disaster 
relief, which guides the budget allocation and positioning of RUAs in order to maximize 


the number of expected survivors. 


The positioning of RUAs has similar challenges to the positioning of facilities 
(see, for example, Shmoys et al. [1997]), particularly the positioning of emergency assets, 


such as fire houses or ambulances, which has evolved since the mid-1960s. The first 


Z 


models proposed were integer linear programming formulations, as noted by Toregas et 
al. [1971], who introduce the “location set covering model” to minimize the number of 
ambulances needed to cover all demand points assuming one type of ambulance and 
unlimited units. Church and ReVelle [1974] propose an alternative approach to maximize 
population coverage subject to limited ambulance availability. Both models assume a 
demand is covered if and only if there exists at least one available unit within a specified 
distance from its location. This assumption is valid if all demands require just one unit to 
respond to them. However, in a relief operation, demands often require response from 
multiple units. In other words, an affected area often requires responses from several 
RUA locations and for different purposes. Alsalloum and Rand [2006] have developed a 
model that looks into identifying the optimal locations of a specified number of 
emergency medical service stations. This model extends the original maximal covering 
location model stated above and considers factors such as people in distress and the 
response time required. These additional factors are taken into account in this thesis, such 
as the need for commodities, transportation capacity, ramp space and different 
transportation means. Hale and Moberg [2005] introduce a selection model for 
emergency resources that can be utilized by logistics managers and supply chain 
continuity teams to determine the appropriate number and locations of storage areas for 


critical emergency equipments and supplies. 


A typical disaster has an unpredictable nature in terms of location, time, and 
magnitude. Therefore it is worth modeling the disaster’s stochastic nature. Liu and Fan 
[2006] use a two-stage stochastic methodology to model the retrofitting of highway 
bridges, where the objective is to minimize the expected system loss. This problem is 
similar to the one addressed in this thesis in terms of the underlying idea: maximizing the 
number of rescued survivors. However, a key difference is that, in their study, Liu and 
Fan look into the strategic resource allocation for critical infrastructure protection and 


hazard prevention, whereas this thesis studies the allocation and expansion of RUAs. 


Many simulation models have been developed to represent how pre-established 
RUAs will respond to disaster relief operations. For example, Marecki et al. [2005] look 


at the interactions between agents and humans during a disaster response using an agent- 


based simulation. This level of detail is catered more toward operational logistics and 


does not determine the positioning of RUAs, which is assumed given data. 


The impetus of this thesis is to develop and solve mathematical optimization 
models to provide guidance in the decision-making process of pre-positioning RUAs for 


disaster relief. 


We claim the selection of optimal position for RUAs and maximization of rescued 
survivors can be guided with a two-stage stochastic optimization model. The contribution 
of this thesis is two-fold: First, the models we have developed integrate the interplays 
among multiple entities: we differentiate between potential survivors in need of 
emergency evacuation and homeless people in need of commodities, and consider their 
spatial location (albeit aggregated by proximity); we account for different transportation 
means and their transportation logistics, including operating time, routing, capacity (for 
workers, commodities and evacuees), and landing needs (if applicable); and we represent 
the location of health personnel, warehouses and ramp spaces, and allow for their budget- 
constrained expansion. We are not aware of any model that explicitly incorporates all of 
these aspects. Second, we account for uncertainty in the problem data by separating 
strategic decisions (to be made now) from operating decisions (which are scenario 
dependent), and demonstrate the benefits of the two-stage approach versus planning for a 
single scenario. While there are multiple models in the area of two-stage stochastic 
programming (some described above), none of them specifically address the problem of 


strategic pre-positioning of RUAs. 


The remainder of this thesis is organized as follows: In Chapter II, we develop 
and describe the mathematical formulation of the pre-positioning model. Chapter II 
describes our test cases and assumptions. Results for these cases are reported in Chapter 


IV. Chapter V provides conclusions and identifies areas for future research. 


I. A MODEL FOR OPTIMIZED PRE-POSITIONING OF RELIEF 
UNITS AND ASSETS 


A. OVERVIEW 
This chapter introduces the Pre-positioning Optimization Model (POM). First, it 
is important to develop a deterministic model which considers the following fundamental 


data: 


e A set of affected areas (AAs) and relief locations (RLs); 


e Potential survivors (PSs) in each AA; 

e Commodities needed in each AA; 

e Workers required to handle commodities in each AA; 

e Available health personnel (HP) and warehouses at each RL; 


e Available ramp space (RS) at each AA; 


e Available transportation means (TM), with associated capacity for 
survivors, commodities and relief workers, time to travel between the 


AAs and the RLs, available operating hours and operating range; and, 


e Allocated budget for pre-positioning of additional HP, warehouses, RS 


and post-disaster engagement of transportation means. 


AAs represent areas hit by the disaster and RLs represent RUA locations, which 
are usually located away from the AAs. In this thesis, it is assumed that PSs are people 
who survive the disaster and need medical evacuation from an AA. For simplicity, all PSs 


are assumed to have the same priority to be rescued. 


PSs are picked up only by TM configured for “special” mission, which depicts a 
TM that can transport injured people but is not configured to deliver commodities. As 
opposed to these, TM configured for “general” mission deliver only commodities and 
relief workers. Each TM’s travel time depends on the distance between the AA and the 
RL, the TM’s speed and the severity of the scenario, which may delay or even impede 


traveling. Land-based TM may not be able to travel to every RL due to damaged or 
5 


impassable roads. Also, some TM are not allowed to use specific RLs, for example, if the 
TM is an aircraft and the RL does not have an adequate runway or a helipad. Each TM is 
capable of performing multiple trips, but it is restricted by its available operating hours. 
Certain TM, such as the CH-53, require ramp space in order to deliver commodities to an 
AA. Land-based TMs, such as_ High-Mobility Multipurpose Wheeled Vehicle 
(HMMWY) (U.S. Army, 2006) do not require ramp space as they are very mobile and are 
assumed to be able to unload their commodities directly to the homeless people in the 


AA. 


It is assumed that every homeless person in an AA is in need of commodities. Not 
meeting this need will cause a deterioration in health, resulting in an increase in the 


number of deaths. In order to capture this feature, we have introduced a penalty function. 


Unloading and organizing the commodities in AAs requires relief workers’ 
intervention. Therefore, commodities are accompanied by relief workers and they share 


capacity on “general” mission TM according to a linear relationship. 


Total cost associated with the expansion of medical facilities and associated HP, 


warehouses, RS, and additional TM is limited by a total allocated budget. 


Stochastic programming provides the platform for modeling optimization 
problems under uncertainty. Our problem is modeled via two-stage stochastic 
optimization. The decision variables of this model are split into two sets. The first set of 
variables (first-stage) is decided before the disaster strikes. The second set of variables 
(second-stage) is based on the concept of recourse: the ability to take corrective actions 
after an uncertain event has taken place. Under the paradigm of a recourse model, we 
have to make some decisions: expansion of warehouses, medical facilities and RS. Then, 
we maximize the expected number of rescued survivors based on the consequences of 
those decisions and the possible values of the random, second-stage variables. These are 
dependent on the scenario which characterizes uncertainty: number of PSs at each AA, 
the demand of commodities which is driven by number of homeless people; and the 


accessibility to the AAs. 


In addition to the above strategic decisions, the model output provides details on 


expenditure of budget, number of rescued survivors, amount of commodity delivered, 
6 


number of relief workers needed, number of TM required and the TM’s approximate 


number of trips between the RLs and the AAs. 


B. PRE-POSITIONING OPTIMIZATION MODEL 

Pre-positioning optimization model (POM) is a two-stage, linear, mixed-integer 
stochastic model. We assume POM looks into pre-positioning decisions in preparation for 
the surge (e.g., within three days after a disaster strikes) of relief operations. Additional 


modeling assumptions are discussed in Chapter III. The POM formulation follows: 


Indices and Index Sets: 
A, Set of affected areas; a € A= {1, 2, 3, ...} 


L, Set of starting and drop off relief locations /; 
VeL= 2.3, es) 


T, Set of transportation means; 


t eT ={CH-53, MV-22, HMMWV} 

i cT, Subset of transportation means that can depart from or drop 
off at relief location / 

TRer ' Subset of transportation means that requires ramp space 


Q, Set of scenarios; @ € O ={1, 2, 3, ae 


Parameters and units: 


Scenario-dependent Data 


dcom®?, 


dsur?, 


Survivor Data 


hpps, 


ichpps ip 


maxehpps it 


vechpps ic 


Ramp Space Data 


icrg, 


MAaXerg , 


Tonnage of commodities needed in affected area a under 


scenario @ [tons] 


Number of potential survivors in affected area a under 


scenario @ [survivors] 


Time taken for transportation means f to travel from relief 


location / to affected area a under scenario @ [hours / trip] 


Relief workers per ton to handle commodities at affected 


area a under scenario @ [workers / ton] 


Probability of scenario @ occurring 


Number of survivors that a health personnel can handle 


[survivors / health personnel] 

Initial number of health personnel of relief location / 
[health personnel] 

Maximum expansion for health personnel at relief location 
/ [health personnel] 

Variable expansion cost for health personnel at relief 


location / [$ / health personnel] 


Initial ramp space capacity at affected area a [tons] 


Maximum expansion for ramp space at affected area a 


[tons for three days] 
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VeClg , 


Commodity Data 





icc), 


maxeCy, 


veccy, 


Variable expansion cost for ramp space at affected area a 


[$ / ton] 


Initial capacity for commodities at relief location / [tons] 


Maximum expansion for commodities at relief location / 
[tons] 
Variable expansion cost for commodities at relief location / 


[$ / ton] 


Transportation Means Data 


ict, 


maxet,, 


vect;, 


comcap;, 


WCaPp; 5 


scap;, 


Initial number of units of transportation means ¢ [number of 


units] 


Maximum expansion for transportation means ¢ [units for 


three days] 


Variable expansion cost for additional unit of transportation 


means ¢ [$ / unit] 


Commodity capacity of “general” transportation means ¢ (if 


loaded with cargo only) [tons / transportation means x trip] 


Relief worker capacity of “general” transportation means ¢ 


(if loaded with relief workers only) 
[workers / transportation means x trip] 


Survivor capacity of “special” transportation means ¢ 


[survivors / transportation means x trip] 


Daily available hours of transportation means ¢ 


[hours / aircraft x day] 


Operating range of transportation means ¢ [hours] 
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qe, 


Miscellaneous Data 





Total budget allocated [$] 


Penalty by unmet commodities (i.e., gc homeless people 
are assumed to perish per ton of unmet commodities) 


[survivors / ton] 


Derived Sets and Data: 


yO 


L°eTE, 


AeA 


\ hea eaed 2 


Vivendi 


Subset of relief locations where survivors could be dropped 


off; 


derived as € Llichpps, >0 or maxehpps, > 0| 


Subset of relief locations from where commodities could be 


supplied; 
derived as € Lliecy >0 or maxecy > 0 


Subset of affected areas where ramp space exits or may 


exist; 

derived as {a E Alicrg > 0 or maxerg > 0} 

Subset of transportation means used for “general” missions; 
derived as \t eT |scap, = 0| 

Subset of transportation means used for “special” missions; 


derived as \t eT |scap, > 0| 
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KReT RE eASL, Subset of four-tuples (¢, J, a, /') where it is feasible for 


transportation means ¢ to travel from / to a and then to /'; 


derived as 


{(Wal')eTxLxAxL yg tty SM te Tr Ane 


K° CK, Subset of four-tuples (¢,/,a,/') e.la,l')eK,teTO el 


where it is feasible for “general” transportation means ¢ to 


travel from / to a and then to/’. 


K° CK, Subset of four-tuples (t,/,a,/') (.1a,!') eK,teT® les 


where it is feasible for “special” transportation means ¢ to 


travel from / to a and then to/'. 
ics) Initial capacity for survivors at relief location / [survivors]; 
calculated as ics; = ichpps; x hpps 


maxes Maximum expansion for survivors at relief location / 


[survivors]; 
calculated as maxes; = maxehpps; x hpps 


vecs| Variable expansion cost for survivors at relief location / 
[S]; 


calculated as vecs; = vechpps;/ hpps 
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Decision variables and units: 


Commodity Decision Variables 


COMD® », 


TCOMD? 


ta? 


UGE. 


HC): 


Commodities delivered by transportation means f¢ traveling 
from / to a and then /7' under scenario @ [tons] 


Total commodities delivered by transportation means f¢ to 


affected area a under scenario @ [tons] 


Unmet commodities at affected area a under scenario @ 
[tons] 
Expansion needed for commodities at drop off relief 


location / [tons] 


Survivor Decision Variables 





NSURR® »», 


TNSURR®, 


ES}, 


Number of survivors rescued by transportation means ¢ 


traveling from / to a and then /' under scenario @ 


[survivors] 


Total number of survivors rescued by transportation means 


t at affected area a under scenario @ [survivors] 


Expansion needed for survivors at drop off relief location / 


[survivors] 


Ramp Space Decision Variables 


ERa, 


Expansion needed for ramp space at affected area a [tons] 
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Transportation Means Decision Variables 





ETM, Additional transportation means ft needed under scenario @ 


[number of units] 


NT. RIP. I! Number of trips from / to a and then to /’ by transportation 


means ¢ under scenario @ [trips] 


TWD, Number of relief workers carried by transportation means ¢ 


to affected area a under scenario @ [workers] 
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Formulation: 


(1) max) p? >» > TNSURR# - ¥gcUC? 
Q@ a 


teTS 4 





Subject to: 

Budget Constraint 
- Total budget: 

(2) De vecs, E'S; + Ds: vecc, EC, + > vecly ERg ee ETM? <b 
leLS leL© ac AR 


Commodity Constraints 





- Maximum expansion: 
(3) EC, <maxec, Vile L& 


- Maximum supply: 
(4) »y COMD® ,, Sicc, + EC, VWieL°,a 





- Route capacity: 


(5) COMD?, , <comcap,NTRIP? |, 


- Total tonnage: 





(6) TCOMDZ= = Y COMD?., VteT?,a,a 
(2 [*/,a,l'jeKe 
- Meet demand: 
(7) > TCOMD® +UC® =dcom? Va,o 
teTG 
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Vt,l,a,l' |(t,,a,l’)eK?,@ 


Transportation Means Constraints 


- Maximum extra units: 
(8) ETM? <maxet, 


- Operating hours: 








(9) (ta +a) NTRIP Ry Shy (ict; + ETM”) 
(/.4,1' (t,l,a,/")eK 
- Balance at RLs : 

(10) > NTRIPP n= NT RIPO 
(a,l)\(t,L.a,l)eK (',a)\(t,La,l")eK 

Survivors Constraints 
- Maximum expansion: 

(11) ES) <maxes, 
- Maximum capacity: 

(12) > NSURR® j, Sics), + ES), 


(1.4) 





(‘.l.a,l')exS 


- Route capacity: 
(13) NSURR?,, <scap,NTRIP® , 
- Total survivors: 
(14) ZTNSURR@ = » NSURR® y 
(u(tLa,"}eKS 





- Maximum survivors: 
(15) >; TNSURRQ < dsur® 


teTS 
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Vt,o@ 


Vt,o@ 


VLteT!,@ 


VleE LS 


Viel ,a 


Vt,l,a,l'|(t,l,a,l') €K*,@ 


VteT*,a,o 


Va, 


Relief Workers Constraints 





- Relief workers follow commodities: 
(16) >) TWD2>wpc® SY) TCOMDE 


teT© teT@ 


- Joint capacity on TM: 
(17) weap,TCOMD? + comcap, TWD? < 
wcap,comcap, > NTRIP? 1 
(LU)[tLa,l"\eKe 





Ramp Space Constraints 

- Maximum expansion: 
(18) ER, <maxer, 

- Maximum capacity: 
(19) }) TCOMDS <icry + ER, 


teTR 


Non - negative Variables : 








(20) COMD? ,, Vt,l,a,l'|(t,l,a,l)eK?%,@ 
TCOMD®, VteT?,a,o 
UC?, Va,o 
EC, Vie LE 
NSURR® ,,, Vt,/,a,!'|(t,l,a,l') €K*,@ 
TNSURR®, VteT*,a,@ 
ES), Vie LS 
ER,, Vae AR 
TWD2, VteT?’,a,o 

Non - negative Integer Variables: 

(21) NTRIP?), Vilal'|\(t,Lal)eK,o 
ETM, VieT/,o 
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Va,o 


Va,t|teT?,@ 


Vae AR 


Vae AR @ 


C. MODEL DESCRIPTION 

1. Objective Function (1) 

POM seeks to maximize the total expected number of rescued survivors, by all 
TM and from all AAs. However, a penalty is applied if unmet commodities at the AAs 
occur. 

2. Budget Constraint (2) 

This constraint ensures the allocated budget is not exceeded. It consists of first- 
stage decisions for expansion of HP, warehouses and RS, and second-stage decisions for 
the use of extra TM from the available fleet. 

3: Commodity Constraints (3)-(7) 

Constraint (3) ensures the allowable warehouse expansion for commodities does 


not exceed the given maximum expansion. 


Constraint (4) limits the commodities that can be delivered from the eligible 


warehouses to the final capacity of the warehouses (after expansion, if any). 


Constraint (5) ensures the commodities carried by each TM (configured for 
“general” mission) traveling on a given route do not exceed the TM’s capacity, which 


depends on the trips the TM makes on that route. 


Constraints (6) and (7) ensure the total amount of delivered commodities does not 
exceed demand, and account for unmet demand, if any. The latter is penalized in the 
objective as described above. 

4. Transportation Means Constraints (8)—(10) 

Constraint (8) ensures additional TM used do not exceed the maximum available 


of each type. 


Constraint (9) ensures the total travel time per TM does not exceed the available 


operating hours of the total (initial plus additional) units of that TM type. 


Constraint (10) is a flow-balance constraint for TM in and out each RL. 
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5. Survivors Constraints (11)—(15) 
Constraint (11) ensures the allowable HP expansion does not exceed a given 
maximum. See the derivation of survivor’s data based on HP data in the “Derived Sets 


and Data paragraph. 


Constraint (12) limits the number of survivors that can be treated in the eligible 
medical facilities (i.e., those located in the respective RLs) based on initial plus expanded 


HP. 


Constraint (13) ensures the survivors carried by each TM, configured for special 
mission, traveling on a given route do not exceed the TM’s capacity, which depends on 


the trips the TM makes on that route. 


Constraints (14) and (15) ensure the number of total rescued survivors does not 
exceed PSs demand at each AA. 

6. Relief Workers Constraints (16) and (17) 

Constraint (16) ensures relief workers arrive at the AAs at a given rate based on 


the amount of commodities supplied to the AA. 


Constraint (17) depicts total capacity on a “general” mission TM as a linear 
relationship between relief workers and commodities. 

Ts Ramp Space Constraints (18) and (19) 

Constraint (18) ensures the allowable RS expansion does not exceed its maximum 


expansion. This constraint is only applicable for RS in AAs where RS exists or may exist. 


Constraint (19) ensures the total tonnage of delivered commodities does not 
exceed the initial RS and the expanded RS, if any, in each AA. 

8. Domains for Decision Variables (20) and (21) 

These constraints define the appropriate domains for the non-negative and non- 


negative integer variables. 
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lil. TEST CASE DESCRIPTION AND ASSUMPTIONS 


This thesis evaluates two notional test cases, which correspond with relief 
operations in an earthquake and a hurricane disaster, respectively. Within each test case, 


two scenarios describe uncertainty. 


A. COMMON DATA AND ASSUMPTIONS 

1. Operation Period 

A typical disaster relief operation will operate throughout days and nights for as 
long as needed until the devastated area returns to normalcy. In this model, we assume 
the only restriction from round-the-clock operation is the respective TM’s daily operating 
hours. In an actual relief operation, the number of working hours maybe hindered by poor 
relief operation planning, long decision cycle from policy makers, aftershocks of an 
earthquake, other logistic delays and social unrest. Typically in a disaster relief operation, 
the first few days after a disaster struck will be the most chaotic. Our scenarios look at a 
hypothetical surge period consisting of the first three days after a disaster. 

2 Setup of Affected Areas and Relief Locations 

We consider a fictitious area which is divided into two AAs—a; and a;— and 
four RLs—1;, />, /; and /,. See Figure 1. The details of each RL and AA will be described 


in the following sections. 


Disaster hit region 





| Relief location region 





Figure 1. Affected areas and relief locations 
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35 Affected Areas 
Three attributes in a typical AA are discussed in this thesis: PS, demand for 
commodities (and associated relief workers) and RS. An overview of our test cases 
follows: 
a. Potential Survivors 
PSs are survivors in need of medical evacuation. For simplicity, we 
suppose every PS in an AA has equal priority to be rescued. For example, a PS who has a 
severe head injury and a PS who has broken his leg shares equal priority to be evacuated. 
We assume the number of PSs in each test case based on the actual fatalities of a disaster 


of that type [Wikipedia, 2005] and [Wikipedia, 2006]. See Table 1. 





Affected| Earthquake test case | Hurricane test case 
area_| Scenario 1 | Scenario 2 | Scenario | | Scenario 2 
ay 70,000 20,000 12,500 20,000 
a2 30,000 80,000 12,500 5,000 
































Table 1. _—_—- Potential survivors by test case and scenario [survivors] 


b. Demand for Commodities and Relief Workers 

We assume the amount of commodities needed is driven by the number of 
homeless people who are not in need of emergency medical evacuation but require food, 
medicines, blankets, etc. Not meeting this demand will result in an aggravated condition 
or even death. See Table 2. In this thesis, we have assumed a penalty factor: gc = 10 


causalities per ton of unmet commodities. 


Relief workers are required to unload and organize the commodities once 
a TM reaches an AA. They are only allowed to travel on TM configured for “general” 


mission, sharing the TM capacity with commodities. 





Affected) Earthquake test case | Hurricane test case 
area | Scenario | | Scenario 2 | Scenario | | Scenario 2 
ay 21,000 6,000 5,000 8,000 
a? 9,000 24,000 5,000 2,000 
































Table 2. Demand for commodities by test case and scenario [tons] 
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c Ramp Space 
RS is required by all air TM. For example, a CH-53 needs RS in order to 


land and deliver commodities to an AA. See Table 3. 


Initial capacity for | Max expansion for | Variable expansion cost 
Affected area | ramp space [tons] ramp space [tons] for ramp space [$/ton] 








ay 2,000 10,000 2,000 
a2 1,000 20,000 1,000 








Table 3. Ramp space attributes for all test cases 


We have assumed some differences in RS at the AAs: 


° Initial RS capacity at a; is twice that of az. 

° Maximum RS expansion for a; is twice that of a). 

° Variable RS expansion cost of a; is twice that of ap. 
4. Relief Locations 


A RL can be an air base, a military headquarters, a warehouse or a medical 
facility with all its components; aircrafts, vehicles, commodities, and doctors. Table 4 
shows the data used for RLs. A discussion follows: 

a. Air Bases and Military Headquarters 

This is the start and drop-off location for any TM. TM can be land-, air-, 
or sea-based. An air base or a military headquarters may have all or a combination of 
TM. For example, during the relief operation, an air base can also be used as a starting 
location for trucks; an open area within a military base can be converted into a landing 
area for helicopters, thus catering to land helicopters. Its medical facilities can also be 
used to treat PSs. In our test cases, we assume all TM are able to start and drop-off at 
every RL, but our model easily accommodates restrictions that may be applicable to other 


Cases. 
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b. Warehouses 

Commodities like food, water, daily human hygiene and medication for 
homeless people are aggregated in this model. These commodities are kept in warehouses 
and cater to the AAs. They can only be transported only by TM configured for general 
mission. 

c. Medical Facilities 

In this model, we have aggregated doctors, nurses, medical facility’s 
administrative personnel, medication requirements, surgery room requirements and 


emergency room requirements as HP. 


Every PS is assumed to be healthy once they are evacuated from AA and 
transported to a medical facility. We assumed each HP can handle up to five PSs. 

d. Differences between Relief Locations 

We assumed four RLs in our test cases. Two of them are located in urban 


areas and two are in rural areas, with the following differences: 
e /; and /> (urban areas): 


- Variable warehouse expansion cost in /; and /, is twice that 


of /3 and 14. 
e /; and /, (rural areas): 


— Initial capacity for both medical facilities and warehouses 


is twice that of /; and />. 


— The maximum expansion for both medical facilities and 


warehouses is twice that of /; and Jp. 


- Medical facility variable expansion cost is twice that of /; 


and J. 
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Earthquake test case Hurricane test case 


















































Relief location Relief location 

Attribute i Lb 1; ly l L 13 lg 
Initial capacity for 
warehouses [tons] 2,000 | 3,000 | 4,000 | 6,000 | 1,000 | 1,500 | 2,000 | 3,000 
Maximum expansion for 
warehouses [tons] 4,000 | 4,000 | 8,000 | 8,000 | 4,000 | 4,000 | 8,000 | 8,000 
Variable expansion cost for 
warehouses [$/ton] 690 | 690 | 345 | 345 | 1,380 | 1,380} 690 | 690 
Initial capacity for medical 
facilities [HP] 300 | 300 | 600 | 600 | 150 | 150 | 300 | 300 
Max expansion for medical 
facilities [HP] 1,000 | 1,000 | 2,000 | 2,000 | 1,000 | 1,000 | 2,000 | 2,000 
Variable expansion cost for 
medical facilities [$/HP] 1,250 | 1,250 | 2,500 | 2,500 | 2,500 | 2,500 | 5,000 | 5,000 





Table 4. Attributes for the relief locations 


=f Transportation Means 

In our test cases, we consider three types of TM. We use two types of 
aircraft (CH-53 and MV-22) and one type of truck (HMMWV). CH-53S and MV-22S are 
configured for “special” mission and CH-53G and MV-22G are configured for “general” 
mission. They are able to handle tasks like rescuing PSs and delivering commodities, 
respectively. HMMWYV can only deliver commodities. See Table 5. 

a. Daily Operating Hours for Transportation Means 

We limit the operation for all TM to 20 hours per day. The remaining four 
hours are assumed for gas refill and change of shift for the flight and land operators. We 
model the three-day surge period by setting total operating hours per TM to 60 hours. We 
also assume that TM are ready to use and no maintenance issues or breakdowns occur. 

b. Travel Time of Transportation Means 

Travel time is calculated based on the TM’s speed and the distance from 
each RL to each AA. For example, a MV-22 is 1.6 times faster, and a HMMWV is five 
times slower than a CH-53. See Tables 6 and 7. /; and /; are assumed to be closer to all 


AAs. 
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Transportation Means 
Attribute CH-53S_ | CH-53G | MV-22S | MV-22G | HMMWV 

Availability [# of units] 20 20 15 2 100 
Maximum expansion [# of units] 20 20 20 20 200 
Variable expansion cost [$ / ton] 500,000 | 500,000 | 800,000 | 800,000 2,000 
Commodity capacity [tons] 0 18 0 18 1 
Survivor Capacity [# of survivors] 24 0 12 0 0 
Worker Capacity [# of workers] 0 55 0 24 8 
Daily available hours [hours] 60 60 60 60 60 
Operating range [hours] 8 8 10 10 14 























Table 5. Attributes for all transportation means 


CH-53 MvV-22 HMMWV 
Affected area Affected area Affected area 


Relief location 


a 
| 60 | 80 | 375 | 50 | 300 | 400 | 
| 80 | 60 | 50 | 375 | 400 | 300 | 


Table 6. ‘Travel time (in minutes) for all transportation means in earthquake test case 














CH-53 MV-22 HMMWV 
(both scenarios) | (both scenarios) 
Relief Affected area | Affected area Affected area Affected area 


location 





Table 7. Travel time (in minutes) a all feangportation means in ee test case 


6. Budget 
“In fiscal year 2004, the average budget for a state emergency management 
agency was $40.8 million, a 23 percent reduction from fiscal year 2003.” [PA, 2005]. We 


use this estimated budget ($40 million) for all of the disaster test cases. 
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Ts Scenario Likelihood 
We assume “Scenario 1” occurs with a probability of 25%, and “Scenario 2” 


occurs with a probability of 75% for both test cases. 


B. DIFFERENCES IN DATA AND ASSUMPTIONS 

The following sections discuss the differences between the two test cases and 
scenarios. 

1. Earthquake Test Case 

In this test case, we suppose the earthquake strikes an urban area. We create two 
hypothetical scenarios and the differences are the number of PSs and the demand for 
commodities in each AA. We also suppose all roads leading to an AA are passable to 


land means in both scenarios. 


For scenario one, the epicenter of the earthquake occurs close to a;, thus a; 
sustains more damage than az. This results in more PSs and higher demand for 


commodities. Scenario two has the opposite situation: more severe damage in ap. 


The total number of PSs is based on the estimated fatalities of the 2005 Pakistan 
earthquake [Wikipedia, 2006]. The demands for commodities are estimated from the 
number of homeless people. 

2. Hurricane Test Case 

In this test case, we suppose a hurricane strikes a rural coastal area. In scenario 
one, we assume both a; and az suffer the same damage. For scenario two, the hurricane 


hits a2 first and then a;. Therefore, az sustains more damage than a). 


The settings for both scenarios are similar to the earthquake scenarios except for 


the following differences: 


e The occurrence of a hurricane is usually more predictable than the 
occurrence of an earthquake and preparations (e.g., evacuation before the 
hurricane strikes) can be made in advance. Therefore, the number of PSs 
and demand for commodities are lower for a hurricane than for an 


earthquake. 
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e Due to flooding, some roads in scenario one (e.g., traveling from /; to a2 


and /> to a;) are impassable to all land means. 


e Severe flooding happened in scenario two; therefore all roads are 


impassable to all land means. 


® The initial capacity of warehouses in the hurricane test case is half of that 


in the earthquake test case. 


e Variable costs of expansion for HP and warehouses in the hurricane test 


case are twice as those in the earthquake test case. 


The number of PSs and demand for commodities information are from Hurricane 


Katrina [Wikipedia, 2005]. 
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IV. RESULTS 


This chapter presents the results to our two-stage stochastic POM for the 
earthquake and the hurricane test cases described in Chapter HI. We also compare the 
stochastic solution with the deterministic wait-and-see bound, to obtain the Expected 
Value of Perfect Information (EVPI) and obtain a heuristic solution which provides the 


Value of the Stochastic Solution (VSS), see Birge and Louveaux [1999]. 


Essentially, the wait-and-see bound represents an upper limit for our model. It is 
obtained by assuming we have perfect information on which disaster is going to happen, 
so we can plan our RUAs accordingly. EVPI is the difference between the stochastic 
solution to POM and the expected value of all individual solutions for each scenario 
under perfect information. EVPI is useful to provide insight on how many more survivors 
could be rescued, on average, if perfect information to plan for the specific disaster were 


available. 


The heuristic solution is obtained by solving a deterministic POM which plans for 
a hypothetical “average” scenario. Values of the first-stage decision variables are then 
fixed in the model to provide different solutions for each actual scenario. VSS is the 
difference between the stochastic solution to POM and the heuristic solution objective 
function. By construction, this solution is sub-optimal, yet still difficult to obtain without 


POM. 


All the computations are executed on an Intel ® Xeon™ CPU, 3.73 GHz 
computer with 3 Gb of RAM running under Microsoft Windows XP operating system. 
The optimization models are coded in General Algebraic Modeling System [Brooke et al. 


1998] and solved by CPLEX 9.0 [ILOG, 2004]. 


In the next sections, our analysis is focused on the strategic decisions regarding 


the allocation of the available budget for expansion. 
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A. EARTHQUAKE TEST CASE 
1. Result Summary 

Figure 2 shows the result for rescued survivors, casualties due to unmet demand 
for commodities, and the total objective function representing a weighed average of the 
two values. The wait-and-see bound yields an optimistic bound (-94,563) on the best 
possible implementable solution (-102,988) provided by our stochastic solution to POM. 
This gives an EVPI of -8,425, or 9% worse than if perfect information were available. On 
the other hand, if our strategic decisions are driven by the “average” scenario, we find a 
heuristic solution (-114,981), which yields a VSS of -11,993, or 12% worse than the 


optimal stochastic solution. 


The difference between the solutions will be explained later with additional 


results. 


Survivors 
Wait-and-see bound Optimal Stochastic Solution Heuristic Solution 
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@ Rescued Survivors O Casualties due to unmet commodities @ Objective value 








Figure 2. Objective Function Breakdown (Earthquake test case). EVPI and VSS 
are 9% and 12% respectively, of the stochastic solution. 
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Figure 3 shows the breakdown of the total budget ($40 million). All budget is 
spent in all scenarios and models. Expenses in warehouse expansion are the least in all 


models, followed by health personnel. 
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Wait-and-see bound 
GB Available @ Used G Health Personnel @ Warehouses O Ramp Space OTransportation Means 
Figure 3. Budget Expenditure Breakdown (Earthquake test case) 
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Table 8 shows the summary of the first-stage decisions on assets that are being 
expanded. 


First Stage Wait-and- Optimal Heuristic 
Decision see bound Stochastic solution solution 
Health Personnel 
(total expanded) 
Health Personnel 
(maximum expandable) 



































(maximum expandable) 
[tons] 























Table 8. Summary of First-stage Decisions (Earthquake test case) 


2. Findings 

The negative values for the objective function shown in Figure 2 are due to the 
heavily penalized unmet commodities: Only 55% and 49% of commodities are delivered 
in the stochastic and heuristic solutions, respectively. Additionally, only 32% and 29% of 
PSs are rescued in the stochastic and heuristic solutions, respectively. Various factors 
such as a high number of PSs, a high demand for commodities and limited numbers of 
TM have contributed to this result. It is observed from Table 8 that all expansions in this 
test case are not under the maximum possible; therefore, budget can be identified as the 


limiting factor. 


We also notice that RS expansion by the heuristic solution is 6% more than in the 
stochastic solution. However, from Figure 3, the heuristic solution has spent 14% less in 
RS than the stochastic solution. This apparent contradiction is caused by the higher 
demand of commodities in a2. We also note that the 14,197 tons of RS expansion by the 


heuristic RS expansion occur only at az. An average scenario is used for the heuristic 
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solution and this drives the solution to be skewed toward optimizing for the second 
scenario. This result shows planning based on an average scenario is less desirable than 


using an optimal stochastic approach. 


B. HURRICANE TEST CASE 

1. Result Summary 

The wait-and-see bound yields an optimistic bound (13,506) on the best possible 
implementable solution (10,082) provided by our stochastic solution to POM. This gives 
EVPI of 3,424, or 25% worse than if perfect information were available. On the other 
hand, using a heuristic solution (-1,348) yields VSS equal to 11,431, or 113% worse than 


the optimal stochastic solution. See Figure 4. 
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Figure 4. Objective Function Breakdown (Hurricane test case). EVPI and VSS are 


25% and 113% respectively, of the stochastic solution. 
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Figure 5 shows the breakdown of the total budget ($40 million) for the hurricane 
test case. All budget is spent in all scenarios and models except by the heuristic model. 
The wait-and-see bound and stochastic solution use most of its budget on the expansion 


of RS whereas the heuristic solution spends most of its budget on TM. 
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Table 9 shows the summary of the first-stage decisions on assets that are being 


expanded. 


First Stage Wait-and- vey Optimal Heuristic 
Decision see bound | “” "* | stochastic solution solution 
Health Personnel 
(total expanded) 
Health Personnel 


(maximum expandable) 


% used 



































(maximum expandable) 
[tons] 























Table 9. Summary of First-stage Decision (Hurricane test case) 


2. Findings 

Figure 4 clearly shows that using a stochastic approach has substantial advantage 
over deterministic solutions. The heuristic solution fails to identify the consequences of 
unmet commodities, although it tends to rescue more survivors than the stochastic 


solution. 
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Vv. CONCLUSIONS AND FUTURE WORK 


We have formulated and solved a two-stage, linear, mixed-integer Pre-positioning 
Optimization Model (POM) to provide guidance in pre-positioning of the Relief Units 
and Assets’ (RUAs). The model takes into account physical limitations, such as the 
capacity of Transportation Means (TM), and other logistic and operational constraints in 
order to maximize the expected rescued survivors and the delivery of required 


commodities over the surge period. 


Two fictitious disaster test cases (earthquake and hurricane) have been tested. To 
account for different levels of damage, we use several scenarios for the number of 
Potential Survivors (PSs) and the demand of commodities in each test case. To add 
further realism, we have also assumed some routes are impassable in some scenarios of 


the hurricane test case. 


In both test cases, the stochastic solution improves the deterministic solution, 
based on the optimal plan for an average scenario, by approximately 11,000 rescued 
survivors. This accounts for almost 11% of the total PSs in the earthquake test case and 


37% in the hurricane test case. 


The available budget is exhausted in both test cases, whereas expansion of 
possible sites is not maximized. This suggests budget is the limiting factor for the deficit 


in the number of rescued survivors and unmet demand in our test cases. 


Future analyses and testing may improve the model and data in order to achieve 


further realism. Some possible areas for research include: 


e Employing real data as estimated by experts in the area of disaster relief 
planning. 
° Testing the robustness of the stochastic solution with respect to each 


individual scenario, and its sensitivity to changes in the data. 


e Incorporating alternative objective functions in POM. Currently, the aim is 
to maximize the expected rescued survivors. Other plausible objectives 


could be minimizing the total budget given a desired level of performance. 
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Incorporating the PSs’ health deteriorating while they wait to be rescued. 
In the current model, we assume there is no health deterioration. The 
model can be further extended to incorporate the survivability curve of a 


PS after a disaster, as well as priorities among PSs. 


POM can be enhanced by incorporating the dynamics of the relief 
operation over time. For example, the delivery of commodities extends 


beyond the surge period at a certain rate until all requirements are met. 


Introducing further survivor classification. There are two classes of 
survivor in this thesis: PSs and homeless people. However, in disaster like 
Hurricane Katrina, survivors can be relocated out of the disaster area. 
Therefore, a survivor can be further classified into another category: 
survivors that need to be relocated to another area, but do not need 


medical attention. 


Changing TM operating patterns. TMs are assumed to travel from a relief 
location (RL) to an AA and back to a RL (not necessarily the same). A 
more reasonable operating pattern is to pick up PSs from different AAs 


(until the TM’s capacity is reached) and then return to a RL. 


Introducing a fixed expansion cost. In order to accurately model the site 
expansion, a fixed expansion cost needs to be introduced. It can be added 
to warehouses, medical facilities and each AA’s ramp space. Similarly, a 


minimum expansion, if any expansion occurs, is sometimes a requirement. 


Providing a graphical user interface to facilitate the input of data and 
display of results. In this way, POM may be used by personnel involved in 


disaster relief planning who are not familiar with optimization techniques. 
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